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Abstract 

We calculate the QCD corrections to the cross section of e + e~ — ► tbH~ and 
its charge-conjugate counterpart within the minimal supersymmetric extension of 
the Standard Model. This process is particularly important if mt < mn + mj, and 
y/s < 2itih, so that t — * bH + and e + e~ — > H + H~ are not allowed kinematically. 
Large logarithmic corrections that arise in the on-mass-shell scheme of quark mass 
renormalization, especially from the tbH~ Yukawa coupling for large values of tan f3, 
are resummed by adopting the modified minimal-subtraction scheme, so that the 
convergence behavior of the perturbative expansion is improved. The inclusion of 
the QCD corrections leads to a significant reduction of the theoretical uncertainties 
due to scheme and scale dependences. 
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1 Introduction 



One of the prime objectives of a future e + e~ linear collider (LC) will be the detailed study 
of spin-zero particles which remain in the physical spectrum after the elementary-particle 
masses have been generated through the Higgs mechanism of electroweak symmetry break- 
ing. Should the world be supersymmetric, then the Higgs sector is more complicated than 
in the Standard Model (SM), which predicts just one neutral CP-even Higgs boson H. 
The Higgs sector of the minimal supersymmetric extension of the SM (MSSM) consists 
of a two-Higgs-doublet model (2HDM) and accommodates five physical Higgs bosons: 
the neutral CP-even h° and H° bosons, the neutral CP-odd A boson, and the charged 
i? ± -boson pair. The 2HDM has six free parameters, which are usually taken to be the 
masses m h o, m H o, m A o, and m H ±, the ratio tan/3 = v 2 /vi of the vacuum expectation 
values of the two Higgs doublets, and the weak mixing angle a that relates the weak and 
mass eigenstates of h° and H°. At the tree level, the MSSM Higgs sector has just two 
free parameters, which are usually taken to be the m^o and tan /3. 

The discovery of the H ± bosons would prove the SM wrong and, at the same time, give 
strong support to the 2HDM and the MSSM. If the H ± bosons have mass nin < ni t — nib, 
they will be mainly produced through the t — > bH + decays of top quarks, which are 
copiously generated singly or in pairs at an e + e" LC [0. On the other hand, if there 
is sufficient center-of-mass (cm.) energy available, <Js > 2m//, then charged-Higgs- 
boson pair production, e + e~ — > H + H~, will be the dominant production mechanism 
|]J. Otherwise, if m# > max(m ( — m&, y/s/2), the H ± bosons can still be produced 
singly. There are various mechanisms of single- charged- Higgs-boson production ||. The 
most important of them are e + e~ — > W + H~, which proceeds through quantum loops 
involving SM |||J and possibly supersymmetric || particles, e + e~ — > t + v t H~ @||, and 
e+e- -> tbH~ §§. 

In the following, we are concerned with the latter process. At the tree level, it proceeds 
through the Feynman diagrams depicted in Fig. [l| It is kinematically allowed if \fs > 
m t +m h +m H . Hence, we are most interested in a situation where m t +m h < y/s/2 < mjj < 
y/8-mt-m b . For yfli = 500 GeV (800 GeV), this implies that 250 £ m H £ 320 GeV 
(400 <^ ran ^ 620 GeV). In such a situation, none of the virtual particles appearing 
in Fig. |l| can be resonating. We note in passing that the absence of resonances is also 
guaranteed if rrit — nib < thh < ni t + m& ||. However, this process is then of minor 
interest because we always have y/s > 2m#, so that e + e~ — > H + H~ will take place. In 
the presence of a resonance, the cross section approximately factorizes. Specifically, we 
have cx(e + e~ — > tbH~) ~ cr(e + e - — > tt)B(t — > bH~) when the virtual top quark gets on 
its mass shell, while we have cx(e + e~ — > tbH~) ~ cr(e + e~ — > H + H~)B(H + — > to) when 
the virtual H + boson gets on its mass shell. 

The purpose of this paper is to investigate the dominant quantum corrections to the 
cross section of e + e~ — > tbH~, which arise from quantum chromodynamics (QCD). The 
relevant Feynman diagrams emerge by attaching one gluon line in all possible ways to 
each of the diagrams shown in Fig. [TJ. This leads to 2 — > 3 diagrams with one closed loop 
(see Fig. 0), which yield the virtual corrections, and to 2 —>■ 4 diagrams of the tree- level 
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Figure 1: Tree-level Feynman diagrams pertinent to the process e + e — > tbH . 



type (see Fig. |3[), which give rise to the real corrections. The loop diagrams involve two, 
three, or four virtual particles. 

While the QCD corrections to the cross sections of the related processes e + e~ — > qqQ, 
where q — t,b and $ = H |7||| or $ = h°,H°,A° ||, are available in the literature, the 
corresponding analysis for e + e~ — > tbH~ has been lacking so far. The present paper fills 
this gap. 

The cross section of the SM process e + e~ — > qqH via a virtual photon and its QCD 
corrections can be recovered from our results as a special case, involving only a subclass 
of the Feynman diagrams shown in Figs. 0-0. As a by-product of our analysis, we confirm 
the numerical results for this cross section obtained in Refs. We also perform the 

complete calculation for this process, which involves Feynman diagrams where the H 
boson is radiated off a virtual Z-boson line, and find good agreement with Ref. 0. 

This paper is organized as follows. In Sec. 0, we list a compact Born formula for the 
cross section of e + e~ — > tbH ~ and give details of our analytical calculation of its virtual 
and real QCD corrections. Lengthy expressions are relegated to Appendices |A| and 
where the Born form factors and the parameterization of the four-particle phase space, 
respectively, may be found. In Sec. we present our numerical results. In Sec. f|, we 
conclude with a summary of our analysis. 
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2 Analytical results 



In this section, we list a compact Born formula for the cross section of e + e~ — > tbH~ 
and give details of our analytical calculation of its virtual and real QCD corrections. By 
charge-conjugation invariance, the results for e + e~ — > tbH + are the same. 



2.1 Born cross section 

We start by defining the kinematics. We call the four-momenta of the incoming electron 
and positron k\ and k 2 and those of the outgoing t quark, b quark, and H~ boson pi, p 2 , 
and p 3 , respectively. We neglect the electron mass, but retain the 6-quark mass, so that 
the on-mass-shell (OS) conditions read k\ = k\ = 0, p\ = mf, p\ = m 2 , and p\ = m 2 H . 
The virtual photon and Z boson have four-momentum p = k\ + k 2 = p\ + p 2 + P3, 
and we define s = p 2 . It is convenient to introduce the dimensionless Lorentz scalars 



a i = Pi I x i = %P ' Pil s -, an d Vi = y %i — 4aj (i = 1, 2, 3). In the cm. frame, Xi = 2p®/ ^fs 
and yi = 2\pi\/y/s carry the meaning of scaled energies and absolute three- moment a, 
respectively. By four-momentum conservation, we have x\ + x 2 + 23 = 2. 
The differential Born cross section may be evaluated as 

1 1 

da-Bom = ^ -\T Born \ 2 dPS 3 (p;pi,p 2 ,p 3 ), (1) 

where the first and second factors on the right-hand side stem from the flux and the average 
over the lepton spins, respectively, 7Bom is the transition-matrix element corresponding 
to the Feynman diagrams of Fig. [I], and the summation over the lepton and quark spins is 
implied. We assume that the incoming leptons are unpolarized. Here and in the following, 
we define the Lorentz-invariant n-particle phase-space measure as 

dPSnfopx, ...,Pn) = (2?r)V 4 ) (p - ±p)i n 7^k-o- ( 2 ) 



i=l 



We now discuss the parameterization of the three-particle phase space. We wish to 
express the Born cross section differential with respect to the scaled energies of the final- 
state quarks, x\ and x 2 . For convenience, we work in the cm. frame, define the z axis of 
the coordinate system to point along fci = — &2, arid fix the x axis in an arbitrary way. 
We then have 

4 

rfPS 3 (p;pi,P2,P3) = j^r^dp l dcos9id(f)idp 2 d(f)2, (3) 

where 9\ and <\>\ are the polar and azimuthal angles of pi, respectively, and <fr 2 is the 
azimuthal angle of P2 with respect to the axis pointing along p x measured from the plane 
spanned by k± and p\. Due to the azimuthal symmetry of the problem at hand, the 
integration over <pi is trivial, and we may take p\ to lie in the x-z plane. If we now rotate 
the coordinate system in such a way that p\ points along the z axis and p^ lies in the x-z 
plane, then 6 = Q\ and <fi = 7r — <p 2 define the direction of fej.. We thus have 

g 

dPS 3 (p;p 1 ,p 2 ,p 3 ) = Tf^^dx! dx 2 d cos 9 d<p. (4) 
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Next, we observe that |7Born| 2 can be written as a contraction of two rank- two tensors, 
a leptonic one involving ki and k 2 and a quarkonic one involving pi, p 2 , and p%. The 
leptonic one has the form 

IT = tr MY (v' e - a! elb ) fo-f(v e - a e75 ), (5) 

where v e , v' e , a e , and a' e are generic vector and axial- vector couplings of the electron to 
the photon or Z boson. Performing the integrations over cos 9 and <p, we obtain 

r dfl 4 

J 4^ = 3 {VeV ' e + GeG ' e) W ~ Sgn ' (6) 

The fact that Eq. (|6]) just depends on p dramatically simplifies the remaining phase-space 
integrations, since scalar products of the type ki ■ pj are precluded. The residual scalar 
products are pi ■ pj = (s/2)(zk — a% — ctj), where Zk = 1 + au — Xk, with k — 1, 2, 3 and 

We thus find the doubly-differential Born cross section to be 



1 



dxi dx 2 327r 3 v / 2 
where Gf is Fermi's constant, 



Q%j(x u x 2 ) + Q e Vefjz(xi, x 2 ) + (v 2 + Al) f zz (x 1 ,x 2 )] , (7) 



h - 2s 2 w Q. 



O — —9r s O V — a 



2 / ' e i 111 \ 

m z /s 1 — m z /s 

), and fzz{xi,x 2 ) are form factors listed in Appendix 0. Here, 
= 1 — (? w = 1 — m^/m| is the sine square of the weak mixing angle, mw and m z 
are the masses of the W and Z bosons, respectively, Q e = — 1 is the electric charge of 
the electron, and I e = —1/2 is the third component of weak isospin of its left-handed 
component. 

The boundaries of integration are 

< xi < 1 + a x — (y/a2~ + v^) 2 , 
x 2 < x 2 < x 2 , (9) 

with 



%2 " 2 Zl 



[2-x 1 )(z 1 + a 2 -o 3 ) ±y 1 J\(z 1 ,a 2} a 3 



(10) 



where X(x, y, z) = x 2 + y 2 + z 2 — 2(xy + yz + zrr) is Kallen's function. We perform the 
integrations over x\ and x 2 numerically with the aid of the multi-dimensional Monte Carlo 
integration routine VEGAS JL0|. 
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Figure 2: One-loop Feynman diagrams pertinent to the process e + e — > tfe/J . They are 
classified according to the number of propagators in the loop. 

2.2 Virtual corrections 

We now turn to the virtual QCD corrections, which arise from the one-loop Feynman 
diagrams shown in Fig. [|. Specifically, they include t- and 6-quark self-energy corrections; 
ttj, tiZ, bbj, bbZ, and tbH~ vertex corrections; and ttryH~ and tbZH~ box corrections. 
These corrections suffer both from infrared (IR) and ultraviolet (UV) divergences. We 
regularize the former by endowing the gluon with an infinitesimal mass, m g . In our case, 
this does not spoil gauge invariance, since the non-abelian nature of QCD does not yet 
emerge at next-to-leading order (NLO). This leads to terms logarithmic in m g , which 
combine with similar terms arising from soft-gluon emission, to be discussed below, to 
give a m 9 -independent result. We establish this cancellation analytically. 
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UV divergences only occur in the self-energy and vertex corrections. We extract them 
using dimensional regularization, with D = 4 — 2e space-time dimensions and 't Hooft 
mass scale \i. They are removed by renormalization. Specifically, we need to renormalize 
the quark masses and wave functions appearing in Tboto- Notice that the quark masses 
enter not only through the quark propagators, but also through the tbH~ coupling. To 



this end, we substitute m q 



m q +5m„ and ip„ — > ipn 



6Z$ 



t, b), where 



bare quantities are denoted by the superscript 0. In the OS scheme, the renormalization 
constants read 



5mq = ~^° F 



5Z$ = -^C F 

47T 



3A + 31n-^ + 4 + 0(e) 



A + ln-^- + 21n-f + 4 + e>(e) 



m 



"q ■■"q 

with C F = (Nl - 1)/(2JV C ) = 4/3 for N c = 3 quark colors and 

A = - - 7B + ln(47r), 
e 



ID 



(12) 



(13) 



where je is the Euler-Mascheroni constant. Notice that Zf is also IR divergent. The 
expression for 5m q in the modified minimal-subtraction (MS) scheme [|ll|] emerges from 
Eq. ( |TTD by retaining only the first term contained within the square brackets. 
The virtual QCD corrections may be evaluated as 



dcr virt — rfCTBorn^virt^g) 



with 



J virt 



Born 



- Re ^ T Boni 



=t,i, 



dm. 



1 

— ( 

2 



+ ^ 



^Born + lit 



loop 



(14) 



(15) 



where T\ Q o-p ls the transition-matrix element corresponding to the Feynman diagrams of 
Fig [2| Notice that the quark masses that appear in the squares of the quark spinors and 
in the boundaries of the phase-space integration correspond to renormalized ones from 
the outset. As mentioned above, 8 v i r t{m g ) is UV finite, but IR divergent. Notice that 
Eq. (Q), which refers to the physical case D = 4, can still be used at the one- loop level, 
since the quarkonic tensor is by itself UV finite upon renormalization. 

We generate 7i oop and reduce it to standard one-loop scalar integrals in two indepen- 
dent ways: one is based on the combination of the program packages FeynArts [1^] and 



FormCalc ||T3| and the other one on custom-made routines written in the program language 
FORM |l4j . We then evaluate the standard one- loop scalar integrals, the IR-divergent ones 
analytically using the results of Ref. and the IR-finite ones numerically with the help 
of the program package LoopTools 



13 1 . Our analytic result for S v 



[rrir. 



is too lengthy 



to be presented here. 
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Figure 3: Tree-level Feynman diagrams pertinent to the process e + e — > tbgH . 



2.3 Real corrections 

We now proceed to the real QCD corrections, which arise from the 2 — > 4 tree-level 
Feynman diagrams shown in Fig. |^. We denote the gluon four-momentum by q. As 
mentioned above, the emission of soft gluons generates IR divergences. For consistency 
with the evaluation of the virtual QCD corrections, the latter must also be regularized 
by the gluon mass m g . It is convenient to work in the cm. frame and to introduce 
an unphysical gluon-energy cutoff, E CVLt , with m g <C E cut <C E max , where -E max is the 
maximum gluon energy allowed by kinematics, so as to separate the gluon phase space 
into soft and hard regions, defined by m g < q° < E cut and E cut < q° < -E max , respectively. 
This has two technical advantages, since soft gluons with infinitesimal mass m g do not 
affect the kinematics of the underlying process, while hard gluons with zero mass do 
not produce IR divergences. On the one hand, the soft-gluon bremsstrahlung may be 
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treated analytically in the eikonal approximation, which is independent of the underlying 
process and results in a multiplicative correction to the Born result. On the other hand, 
m g may be safely neglected in the treatment of the hard-gluon bremsstrahlung, which 
facilitates the phase-space integration. In turn, the soft- and hard-gluon contributions 
both depend on E cut , while their combined contribution is, of course, independent of E cut , 
which we checked numerically. As mentioned above, the m g dependence of the soft-gluon 
contribution analytically cancels against the one of the virtual QCD corrections. 
The soft-gluon contribution is given by 



(16) 



with 



#soft(^g, -Ecut) 




Pi 



P2 



\g\<E cut q" \pi -q Vi-q, 



(17) 



where q° = Jq 2 + m 2 is the gluon energy. The integration in Eq. (]17|) can be performed 



analytically as described in Ref. [16], the result being 



'a z 3 — ai — a 2 , a 2 ai _\ l 4:E 2 ut 



71 



2 a 2 ai — 0,2 



In 



- 1 In 



(1-2 



+ oi 



z 3 ~ &l ~ 0-2 



a 2 a\ 



Li 2 1 



- a 2 

X2 + V2 



Li 2 fi_ a ^i±ll') + Li 2 (l 



v 



a- 



v 



- Li 2 1 



— In — 



V\ ^2 X 2 

1 In ■ 



X2 - V2 



V2 



+ 1 ln 2 ^ 
4 V + yi 



In 



2 X 2 ~ V2 
%2 + V2 t 



a \y x x x +y! y 2 x 2 + y 2 , 
where Li 2 (x) = — Jq dy ln(l — y)/y is the Spence function, 



(18) 



a 



2ai 



z 3 - ai - a 2 + J X(z 3 , oi, a 2 



a*(Xi — a 2 

v = 2 . 

axi — x 2 

Notice that Eq. ([18]) is invariant under the interchange of the indices 1 and 2. 
The hard-gluon contribution may be evaluated by integrating 



(19) 



^Chard(-E'cut - 



1 1 

2^ 4 



^eal| 2 %° - E cnt )dPS^p;Pl,P2,P3,q), 



(20) 



where % ea i is the transition-matrix element corresponding to the Feynman diagrams of 
Fig. |3], over the full four-particle phase space, imposing the condition q° > E cnt . We use 
the parameterization of the four-particle phase space presented in Appendix [B] It involves 
five nontrivial integrations, which we perform numerically using the Monte Carlo routine 



VEGAS HlO| . Our formula for |7^ ea i| is too lengthy to be listed here. 
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We performed several checks for our implementation of the four-particle phase-space 
integration. We numerically verified the analytical formula for the total cross section of 
e + e~ — ► qqg* — > qqQQ, where q and Q represent massless and massive quarks, respec- 



tively, and g* denotes a virtual gluon, given in Eq. (2) of Ref. [17]. In this case, IR 
singularities do not appear in intermediate steps, so that no separation into soft-gluon 
and hard-gluon contributions is required. We also found excellent agreement with a nu- 
merical result for a similar process involving four different quark masses obtained using 
the democratic multi-particle phase-space generator RAMBO ||18|| . 

Our final result for the QCD-corrected differential cross section reads 

d<?QCD = dcr BoTn [l + 5 virt (m g ) + 5 soit (m g , E cut )} + da h3uTd (E cnt ) , (21) 

where da Bo m, S vixt (m g ), 5 solt (m g , E cut ), and da havd (E cut ) are defined in Eqs. (0), (||), ([H]), 
and (]2"0"D, respectively. It is manifestly independent of m g and insensitive to the choice of 
E cnt , as long as m g -C E cnt -C E m3jX , as we verified numerically. We also checked that the 
QCD-corrected total cross section is finite in the limit m b — > 0, in compliance with the 



Kinoshita-Lee-Nauenberg theorem [IS 



The QCD-corrected cross section of the SM process e + e~ — > qqH via a virtual photon 
can be obtained from our results as a special case, involving only a subclass of the Feynman 
diagrams shown in Figs. |l]-[| As a by-product of our analysis, we confirmed the numerical 
results for this cross section obtained in Refs. [0,f|. We also performed the complete 
calculation for this process, which involves Feynman diagrams where the H boson is 
radiated off a virtual Z-boson line, and found good agreement with Ref. In turn, this 
provides a nontrivial check for all parts of our analysis. 



3 Numerical results 



We are now in a position to present our numerical results. We first specify our input 
parameters. We use m w = 80.419 GeV, m z = 91.1882 GeV, m t = 174.3 GeV, m b = 
4.6 GeV, G F = 1.16639 x 10~ 5 GeV~ 2 [pO] , and the present world average af\m z ) 



0.1180 [|2T 



We consistently evaluate aT'\^) and rUq (//) to lowest order (LO) in the 
MS scheme with nj = 6 active quark flavors performing the matching with n/ = 5 QCD 
at scale m t . For the reader's convenience, we collect the relevant formulas here [22j: 



a 



(5) 



(A*) 



(5)/ 

as '{mz, 



+ 



In 



(6) 



a s ' K 



+ 



7T 

(6) 




2 ' 



7T 



m 



1 V 2 
In—, 

mf 



(6) 



m 



(ji) = mt 



1 - 



a 



(6) 



{mt) 



TV 



a 



(5) 



(m b ) 



C h 



71 



«i 6) (/^) 1 






oa (rn t ). 


! 








af\m t ) 


70/^ 5) 




a s [m h )_ 




af\m t ). 



jo/ft 



(6) 



:22) 



(23) 



(24) 



(25) 
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where 



7o = \c r% (26) 

with Ca = N c and Tp = 1/2, are the first coefficients of the Callan-Symanzik beta function 
and the quark-mass anomalous dimension, respectively. For simplicity, we use a common 
renormalization scale \i in af\ii) and rn^(fi). We study the cases ^Js = 500 GeV 
and 800 GeV. As for the MSSM input parameters, we consider the ranges 1 < tan j3 < 
40 « m t /m b and 250 < m H < 320 GeV if y/s = 500 GeV or 400 < m H < 620 GeV if 
= 800 GeV. 

We now discuss the influence of the QCD corrections on the total cross sections of 
e + e~ — > tb~H~ and its charge conjugate counterpart, which we add. We start by selecting 
the renormalization scheme and scale that are most appropriate for the problem under 
consideration. For this purpose, we study the /i dependence of the Born and QCD- 
corrected results in two different renormalization schemes. The first one uses the pole 
masses m t and m& as basic parameters (OS scheme), while the second one uses m t and 
the MS mass Tn^\jj) as basic parameters (mixed scheme). Both schemes employ the MS 
definition of a^(//). We refrain from utilizing m[ (fi), which, in general, significantly 
deviates from m t , as may be seen from Eq. ([24]) . For, if we were to include the weak decays 
of the t and t quarks in our analysis, then, during the propagation of these quarks between 
their production and decay vertices, configurations near their physicalm&ss shells would be 
kinematically favored. As a matter of fact, the experimentally measured invariant masses 
of their decay products are very close to m t . In turn, the phase space of e + e~ — > tb~H~ 
would undergo a significant, yet artificial change of size if it were parameterized in terms 
of mf (//) rather than mt- On the other hand, the use of mjf^/i) is predicated on the 
grounds that it automatically resums large logarithmic corrections that arise if the tbH" 
Yukawa coupling is expressed in terms of mb- A similar feature is familiar from the 
H — > bb decay in the SM ||23|| . This effect is particularly pronounced for large values 



of tan/3 because the tbH~ Yukawa coupling is then approximately proportional to the 
6-quark mass. 

As a typical example, we consider in Fig. |] the case of ^fs = 500 GeV, tan (3 — 40, and 
mil = 260 GeV. We allow /i to vary over two orders of magnitude, from 10 to 1000 GeV. 
In the OS scheme, the Born result is \x independent, while the Q CD-corrected one depends 
on fi via af^jj). In the mixed scheme, the \i dependence enters at LO via m[ 6 ^(/i) and 
at NLO via af\fi) and mf\fi). Obviously, the theoretical uncertainties due to scheme 
and typical scale variations are significantly reduced as we pass from LO to NLO. On 
the one hand, the OS-scheme to mixed-scheme ratio is brought down to the vicinity of 
unity, from 1.46-2.92 to 0.78-1.43, depending on the value of //. On the other hand, the [i 
dependence within the mixed scheme is reduced by a factor of 5, from 0.020 fb to 0.004 fb 
in absolute terms. Furthermore, we observe that, in the OS scheme, the QCD corrections 
lead to a dramatic reduction of the cross section, by 36-67%. As explained above, this 
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0.06 




H(GeV) 



Figure 4: Total cross section of e + e~ — > tbH~ ,tbH + as a function of fi for yfs = 500 GeV, 
tan /3 = 40, and ran — 260 GeV. The dashed and solid curves correspond to the Born and 
Q CD-corrected results, respectively. The upper and lower sets of curves refer to the OS 
and mixed schemes, respectively. 



is because they contain large logarithmic terms of the form a^(fx) In (M 2 /m^), where 
M is a generic mass scale in the ball park of some suitable average of the final-state- 
particle masses, nib, nit, and ran- In the mixed scheme with fi of order M, such terms are 
shifted from the QCD corrections to the Born result, where they are absorbed into the 
running of (p) from jj, = nib to /i = M. This is reflected in Fig. [| by the fact that, 
in the mixed scheme, the QCD corrections are relatively modest, ranging from —39% to 
+30%. Unless otherwise stated, we shall henceforth work in the mixed scheme, which, for 
plausible values of /i, is superior to the OS scheme as far as the convergence properties 
are concerned. 

Let us now turn to the question of how to fix the value of fi in a reasonable way. Scale- 
setting procedures frequently discussed in the literature include the concept of fastest 
apparent convergence (FAC) |24|], the principle of minimal sensitivity (PMS) ||25|| , and 



the proposal by Brodsky, Lepage, and Mackenzie (BLM) [f26| to resum the leading light- 
quark contribution to the renormalization of the strong-coupling constant. The latter 
is not yet applicable to the problem under consideration, which is of LO in the strong- 
coupling constant. The FAC and PMS prescriptions lead us to select the values of /i 
where the Born and QCD-corrected results intersect and where the latter exhibits a local 
extremum, respectively. We observe from Fig. [| that these two fi values approximately 
coincide, at about 60 GeV. Incidentally, in the close vicinity of these two fi values, also 
the QCD-corrected results in the OS and mixed schemes cross over, so that also the 
scheme dependence at NLO vanishes in this neighborhood, at least as for the two schemes 
considered here. Since the /i dependence is logarithmic, a democratic way of combining 
the three scales m^, ra t , and ran is by taking their geometric means, a. = ^fm^hjnu- In 



12 



the present case, this educated guess yields /i ~ 60 GeV, which nicely agrees with the 
triply distinguished point identified above. We checked that this choice works similarly 
well for the case of \/s = 800 GeV, tan (3 = 40, and m H = 410 GeV. We shall henceforth 
employ it, with the understanding that Fig. f| provides us with a useful estimate of the 
theoretical uncertainties due to scheme and typical scale variations, both at LO and NLO. 




M H (GeV) 

Figure 5: Total cross section of e + e" — > tbH~ ,tbH + without (dotted curves) and with 
(solid curves) QCD corrections in the mixed scheme as a function of m# for \fs = 500 GeV 
and various values of tan/3. The middle, lower, and upper sets of curves correspond to 
tan/3 = 2, 6, and 40, respectively. 

Figures || and || refer to \fs = 500 GeV, while Figs. [7| and || refer to \fs = 800 GeV. 
We investigate the m# dependence for various values of tan (3 in Figs. |] and [7| and the 
tan p dependence for typical values of m H in Figs. |] and |[ In each figure, we present 
the Born and QCD-corrected results in the mixed scheme with \i = ^/m^mtmn- For 
comparison, in Figs. || and [8|, we also present the corresponding results in the OS scheme 
with the sam e scale choice. We observe that the total cross sections exhibit minima close 
to tan/3 « \fnhj m b ~ 6, independently of order and scheme. This may be understood 
by observing that the average strength of the tbH~ coupling, which is proportional to 
J ml cot 2 (3 + ml tan 2 /3, is then minimal J27|. Depending on y/s, tan j3, and m#, the QCD 
corrections may be of either sign. By construction, they are generally rather modest in 
the mixed scheme, although they may reach a magnitude of 50% for specific values of 
y/s, tan/3, and mjj, as may be seen from Fig. [?]. On the other hand, in the OS scheme, 



13 




Figure 6: Total cross section of e + e~ — > tbH~ ,tbH + without (dotted curves) and 
with (solid curves) QCD corrections as a function of tan (3 for y/s = 500 GeV and 
ran = 260 GeV. The upper and lower dotted curves refer to the OS and mixed schemes, 
respectively. The two solid curves referring to the OS and mixed schemes lie on top of 
each other. 



the QCD corrections lead to a substantial reduction in cross section at large values of 
tan j3, by up to 50%. As explained above, this may be attributed to large logarithms 
arising from the tb~H~ Yukawa coupling. Finally, we notice that Figs. |] and |8| support the 
observations made in connection with Fig. f|. In fact, owing to our judicious scale choice, 
the Born and QCD-corrected results in the mixed scheme and the QCD-corrected result 
in the OS scheme all approximately coincide, which nicely demonstrates the perturbative 
stability in the mixed scheme and the feeble scheme dependence at NLO. By contrast, 
the perturbative stability in the OS scheme is rather poor at large values of tan f3. 

It is interesting to investigate the relative importance of the contributions due to 
photon and Z-boson exchanges. To this end, we evaluate the photon-induced part of the 
total cross section by putting V e = A e = and compare it with the full result. We find 
that the bulk of the total cross section is due to photon exchange. In fact, for the typical 
values y/s = 500 GeV and m H = 260 GeV, the photon-induced part exhausts 78%, 80%, 
and 82% of the full result if tan/3 = 2, 6, and 40, respectively. 
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4 Summary 

We considered the process e + e~ — > tbH~ and its charge-conjugate counterpart in the 
MSSM, which are among the dominant charged- Higgs-boson production mechanisms at 
a future e + e~ LC if m t + m h < < rriu < y/s — m t — m&, so that t — > 

and e + e~ — > H + H~ are not allowed kinematically. For -y/s = 500 GeV and 800 GeV, 
this corresponds to the mu windows 250 <^ rrin & 320 GeV and 400 ^ mn ^ 620 GeV, 
respectively. We presented a compact Born formula for the cross sections of these processes 
and evaluated their dominant radiative corrections, which arise from QCD. 

We regularized the IR singularities by introducing an infinitesimal gluon mass m g and 
the UV ones by using dimensional regularization. The IR singularities cancelled when the 
virtual and soft real QCD corrections were combined and the UV ones upon renormalizing 
the masses and wave functions of the quarks in the Born transition-matrix element. We 
established these cancellations analytically. We separated the soft-gluon and hard-gluon 
contributions by introducing an unphysical gluon-energy cutoff E cnt and verified that their 
sum is insensitive to the choice of E cut as long as m g «C E cut E max is satisfied. 

We worked in a mixed renormalization scheme, where the strong-coupling constant 
and the 6-quark mass are defined in the MS scheme, while the t-quark mass is defined in 
the OS scheme. In this way, large logarithmic corrections that arise if the tbH~ Yukawa 
coupling is expressed in terms of the 6-quark pole mass are automatically resummed, so 
that the convergence behaviour of the perturbative expansion is improved. On the other 
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hand, the use of the t-quark MS mass does not entail such an improvement, but it rather 
appears somewhat unnatural from the physical point of view. The inclusion of the QCD 
corrections significantly reduces the theoretical uncertainties due to scheme and typical 
scale variations. We found that the QCD corrections to the total cross section may be of 
either sign, depending on the values of y/s, tan/3, and m H , and that they may reach a 
magnitude of up to 50%. 

The e + e~ LC TESLA, which is being developed at DESY, has a design luminosity of 
3.4 x 10 34 cm-V 1 (5.8 x 10 34 cm^s^ 1 ) at y/s = 500 GeV (800 GeV), which corresponds to 



340 ib" 1 (580 fb" 1 ) per year |28j. Thus, a total cross section of typically 0.03 fb (0.07 fb) 
will yield about 10 (40) signal events per year. 

As a by-product of our analysis, we confirmed the numerical results for the QCD- 
corrected total cross section of the SM process e + e~ — > qqH obtained in Refs. 
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A Born form factors 

In this Appendix, we list compact expressions for the form factors / 77 (xi, x 2 ), / 7 ,z(xi, x 2 ), 
and fzz(x x ,x 2 ) appearing in Eq. ([?]). It is possible to combine the propagators of the t 
quark, b quark, and H~ boson with their couplings to the photon and Z boson by defining 
the effective couplings 



where q = t,b and we have identified x t = x x and x b = x 2 . The tbH coupling introduces 
tan f3 dependence through the combinations 



We find 

fzz = 2V t 2 {4a 1 a 2 (2 + 2a x - x x ) + T + [(l + 2 ai )(l + a x + a 2 - a 3 - x x ) - x 2 (l - x x )}} 
+ AV t V b {2a x a 2 (A + 2ai + 2a 2 - 2a 3 - x x - x 2 ) 

+ T + [(l + a 2 - a 3 - xi)(l + ai - a 3 - x 2 ) + (a x + a 2 )(l + ai + a 2 - a 3 ) - aia 2 ]} 
+ 2V b 2 {4a!a 2 (2 + 2a 2 - x 2 ) + T + [(l + 2a 2 )(l + a x + a 2 - a 3 - x 2 ) - x x {l - x 2 )}} 

- 4VtAXl[(l - x 2 )(l + 2ai - x{) — a x + a 2 — a 3 ] 

+ AV t AbT_[(l + a 2 - a 3 - Xi)(l + a x + 2a 2 - a 3 - x 2 ) 
+ (a x — a 2 )(l + a x + a 2 — a 3 ) — 3aia 2 ] 

- 4V b A t T_[(l + a x - a 3 - x 2 )(l + 2a x + a 2 - a 3 - Xi) 

- (ai — a 2 )(l + ai + a 2 — a 3 ) — 3aia 2 ] 

+ 4V b A b T_[(l - xi)(l + 2a 2 - x 2 ) + a x - a 2 - a 3 ] 

- 2A^{4:a x a 2 (2 + 6a x — 3xi) — T + [(l — 6ai)(l + a x + a 2 — a 3 — xi) — x 2 (l — Xi)]} 

- AA t A b {2a x a 2 (A + 6a x + 6a 2 — 6a 3 — 3xi — 3x 2 ) + T + [(l + 2a x + a 2 — a 3 — xi) 
x (1 + a x + 2a 2 - a 3 - x 2 ) - (a x + a 2 )(l + a x + a 2 - a 3 ) + 7aia 2 ]} 

- 2A 2 b {Aa x a 2 {2 + Qa 2 - 3x 2 ) - T + [(l - 6a 2 )(l + a x + a 2 - a 3 - x 2 ) - x x {\ - x 2 )}} 
+ 2V t Hz{Aa x a 2 (l - 2a x + 2a 2 - 2a 3 - x 2 ) + T+[(l + 2a 2 - 2a 3 - Xl ) 

x (1 — ai — a 2 — a 3 - xi - x 2 ) - 2a x (a x + a 3 ) + 2a 2 (l + ai + 2a 2 - 2a 3 )]} 

- 2V b H z {Aa x a 2 (l + 2a x - 2a 2 - 2a 3 - x x ) + T+[(l + 2a x - 2a 3 - x 2 ) 

x (1 — a x — a 2 — a 3 - x x — x 2 ) - 2a 2 (a 2 + a 3 ) + 2ai(l + 2a x + a 2 - 2a 3 )]} 

- 2AtHzT-[(l + 2ai + 2a 2 - 2a 3 - xi)(l + a x — a 2 — a 3 — x x — x 2 ) 




l-Xq 



(A.l) 



T± = a x cot 2 (5 ± a 2 tan 2 /3. 



(A.2) 
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+ 2a 2 (l - 2ai + 2a 2 - 2a 3 )] 

- 2 </ 4 6 7^ z T_[(l + 2ai + 2a 2 - 2a 3 - x 2 )(l - ai + a 2 — a 3 — x\ — x 2 ) 
+ 2ai(l + 2ai - 2a 2 - 2a 3 )] 

- H|[4aia 2 (3 - 4a 3 - 2x t - 2x 2 ) 

+ T + (3 - 4a 3 - 2xi - 2x 2 )(l + a\ + a 2 - a 3 - xi - a; 2 )]. (A. 3) 

The formulas for f 11 {x 1 ,x 2 ) and f 7 z(xi,x 2 ) m &y be obtained from Eq. ( |A.3| ) by adjusting 
the coupling constants. Specifically, one has to substitute 

V q ^Q q , A^O, Hz^H,, (A.4) 

in the first case and 

V q V q > -> Q ? V g > + Q 9 >V ? , V ? A?' -> Q ? ^, A A -> 0, 

V q 1~Lz — >■ + V 9 7^ 7 , — >■ -A.qli.~f, H 2 Z — > 2H~fHz (A. 5) 

in the second one (g, g' —t,b). 

B Four-particle phase space 

In this Appendix, we present the parameterization of the four-particle phase space that we 
use to evaluate the hard-gluon contribution. We generically denote the four-momenta and 
masses of the final-state particles as pi and mi (i = 1, . . . , 4), respectively. Similarly as in 



Sec. we define p = J2t=iPi, s = p 2 , ai = m 2 /s, Xi = 2p-pi/s, and %ji = Jx 2 — ^a^. Due 
to four-momentum conservation, we have J2f=i Xi = 2. We decompose the four-particle 
phase space into three nested two-particle phase spaces as E[ 



dPS 4 (p;pi,p 2 ,p 3 ,p 4 ) = ^^y^ dPS 2 (p; p 12 , p 34 ) dPS 2 (p 12 ;p x ,p 2 ) dPS 2 (p u ; p 3 , p 4 ), (B.6) 



cm. frame, take the z axis of the coordinate system to point along k\ = — fc 2 , and choose 



where pij = pi +pj and s^- = p?-, with = (1, 2), (3, 4). As in Sec. we work in the 



the x axis arbitrarily. We have pi 2 = — P34 and |pi 2 | = (l/2)wA(s, si 2 , s 34 )/s. Using 



dPS 2 (p;pi 2 ,p 34 ) = (jcos6>i 2 d0i2, 



P12I 

(4tt) 2 VT 

rfPS 2 (pij;Pi,Pj) = ( 47r )2| p ..| rf Pi ( B - 7 ) 

where 6 l i 2 and 0i 2 are the polar and azimuthal angles of pi 2 , respectively, and fa is the 
azimuthal angle of pi with respect to the axis pointing along pij measured from the plane 
spanned by k± and p^, we obtain 

g 

dPS 4 (p;pi,p 2 ,p 3 ,p 4 ) = =dsi 2 dsu d cos 9 12 dfa 2 dp® dfa dpi d<fi 3 . (B.8) 

(4vr)yA(s,Si 2 ,s 34 ) 
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Owing to the azimuthal symmetry of the problem under consideration, the integration 
over 12 is trivial, and we may choose px2 to lie in the x-z plane. If we now rotate the 
coordinate system in such a way that p±2 points along the z axis and p± lies in the x-z 
plane, then 9 = 9\ 2 and (f) = tx — <pi define the direction of k\. Introducing zij = Sy/s, we 
thus have 

s 2 

dPS^p; Pi, P2, P3, Pi) — / ==dz-[2 dzu dxi dx 3 d(f> 3 d cos 9 d(f>. (B.9) 



(47r)yA(l, 2:12,2:34) 

As explained in Eq. (H), the integrations over cos 9 and can be exploited to transform the 
leptonic tensor L^ v defined by Eq. (|), which depends on k\ and k 2 , into one depending 
only onp, so as to preclude scalar products of the type ki-pj. The residual scalar products 
read 

s , 



Pi 


■P2 


P3 


■Pi 


Pi 


•P3 


Pi 


Pi 


P2 


•Ps 


P2 


■Pi 



(234 — ^3 — O'i)-, 

[xix 3 — yx yz (sin 9\ sin 9 3 cos (p 3 — cos 9i cos #3 



[xi - Z12 — di + a 2 ) 



Pi ■ P3, 



(x 3 - z 34 -a 3 + a A ) - px ■ p 3 , 
(l-xi-x 3 + a 1 -a 2 + a 3 - a 4 ) + pi ■ p 3 , 



where 9i is the angle enclosed between pij and pi. It is determined by 



cos 9j 



Xi(l + z i:j - z M ) - 2(zij + cii - dj) 



, Zjj , Zkl ) 

with (k,l) = (1,2), (3,4) and (i, j) 7^ (k,l). Furthermore, we have 

Xj 1 X{ ~\~ Zij Zfc[. 

The limits of integration are 

(VdT+ v^) 2 < 2l 2 < (1 - 
(y/d^ + y/d~l) 2 < 2 34 < (1 - 

< 03 < 2tt, 



(B.10) 



[B.IV 



;b.i2) 



a 3 - 
-12J 2 , 



where 



x 



± 



22, 



(1 + Zij - z k i)(zij + a { - dj) ± ^X(l,Zij 1 z k i)X(zij 1 a h a j ) 



(B.13) 

;b.i4) 



For the application in Sec. |2j.3, it is convenient to assign the indices 1, 2, 3, and 4 to 
the t quark, b quark, gluon, and H~ boson, respectively. Then, the hard-gluon condition 
q° > E cut may be implemented by substituting x 3 —>■ max (x 3 ,2E cut / in Eq. ( |B.13|) . 
Furthermore, we have a 3 = throughout this appendix. 
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